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Abstract 

f^ We solve the optimal trade execution problem in the Almgren and Chirss 

Li^ framework with the Value-at-risk / Expected Shortfall based criterion of 

Gatheral and Schied when the underlying unaffected stock price follows a 
displaced diffusion model. The displaced diffusion model can conveniently 
model at the same time situations where either an arithmetic Brownian mo- 
tion (ABM) or a geometric Browinan motion (GBM) type dynamics may 
prevail, thus serving as a bridge between the ABM and GBM frameworks. 
^ We introduce alternative risk criteria and we notice that the optimal trade 

CN execution strategy little depends on the specific risk criterion we adopt. In 

^^ most situations the solution is close to the simple Volume Weighted Average 

^vq Price (VWAP) solution regardless of the specific diffusion dynamics or risk 
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criterion that is chosen, especially on realistic trading horizons of a few days 

^^ or hours. This suggests that more general dynamics need to be considered, 

^^ and possibly more extreme risk criteria, in order to find a relevant impact on 

. . the optimal strategy. 
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1 Introduction 

A fundamental problem for algorithmic traders is how to partition a large trade 
into smaller trades so as to minimize the effect of market impact. In a number of 
seminal papers including [T] and |2], Almgren and Chriss formalized this problem 
by combining expected execution cost and risk. The framework in [1] and [2] has 
linearly increasing execution costs (in the trading rate) whereas the risk criterion 
is given by the variance. The advantage of such a setup is the possibility to find a 
closed-form analytical solution to the optimal execution problem. This solution is 
static and is usually in the class of the Volume Weighted Average Price (VWAP) 
solutions. A disadvantage of this framework is that it needs to assume an extremely 
stylized dynamics for the underlying unaffected stock, namely an arithmetic Brow- 
nian motion (ABM), possibly leading to arbitrarily negative asset values. Another 
consequence of the framework is that adaptive trade schedules formulated on the 
current stock price level are suboptimal. 

Gatheral and Schied [3] solve the same problem under the more realistic as- 
sumption that the unaffected underlying stock follows a geometric Brownian motion 
(GBM). Besides being the cornerstone of modern basic option pricing theory, GBM 
has the clear advantage of implying positive stock prices in all scenarios. However, 
under the GBM assumption with variance as risk criterion, the optimal execution 
problem is not tractable. The problem has been studied numerically in [4j and [5], 
where it is found that under GBM the optimal trading rate increases with the stock 
price and that the cost-risk efficient frontier is almost the same under ABM and 
GBM, implying that the static solution obtained by solving the ABM version is not 
very suboptimal even under the different GBM assumption. 

However, Gatheral and Schied look for an analytic solution rather than a nu- 
merical one. This leads them to resorting to a different risk criterion, namely the 
Value at Risk measure. This is adopted in [3J and the problem is solved analytically. 
An explicit closed form expression for the trading rate is obtained. The authors are 
then able to analyze precisely the optimal strategy under the GBM assumption, and 
can explicitly compare the optimal strategy under GBM with the optimal strategy 
for the risk criterion under ABM. 

In this paper we extend the Gatheral and Schied [3J analytic result to a more 
general dynamics, namely the displaced diffusion (DD) model. This dynamics is 
more general than both ABM and GBM, but comprises behaviours very similar to 
either ABM or GBM for different values of the parameters in the dynamics itslef. 
In particular, it allows to model volatility as a combination of constant absolute 
volatility (as in the ABM) and level-proportional absolute volatilty (as in GBM). 
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This is a step forward in that it can be seen as a model that unifies the Almgren 
and Chriss and Gatheral and Schied frameworks into a single tractable setup. 

We further introduce a risk criterion we call SAE (Squared-Asset Expectation) 
to investigate how a different risk criterion would change the optimal solution. We 
solve the related problem through an approximation related to calculus of variations 
and we notice that the optimal trade execution strategy seems to depend little on 
the specific risk criterion we adopt. In particular, we find that in most situations the 
solution is close to the Volume Weighted Average Price (VWAP) solution regardless 
of the specific diffusion dynamics or risk criterion that is chosen. This holds unless 
extreme and unrealistic assumptions are considered either for the dynamics param- 
eters or for the risk criterion parameters. This leads to suggest that more general 
dynamics need to be considered, and possibly more extreme risk criteria, in order 
to find a relevant difference in the optimal strategy. 

The paper is organized as follows: Section [2] introduces the optimal trade execu- 
tion problems with a few examples of stylized dynamics for the unaffected asset price, 
and with simple market impact and cost assumptions. In addition, risk functions 
are introduced and Value at Risk (VaR) and Expected Shortfall (ES) are discussed 
in detail. Section |3] introduces the optimal solution for the case where the dynamics 
is a displaced diffusion and where the risk function is either VaR or ES. Section |4] 
explains how the DD with VaR/ES criteria escapes Schied's invariance result. Sec- 
tion |5] introduces a different Squared-Asset Expectation (SAE) risk criterion that 
is closer to the variance, and derives approximate solutions of the execution prob- 
lem under that criterion. Section |6] compares the optimal execution solutions under 
different criteria and dynamics, with numerical examples and discussion, whereas 
Section [7] concludes the paper. 

2 The trade execution problem under different 
dynamics 

Following pp, [2] and especially [3] we assume that the trader's number of shares 
follows an absolutely continuous trajectory 

t \-^ x(t). 

Given this trading path, the price at which the transaction occurs is given by 

St = St + T] x{t) + 7(x(t) - Xq) 
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where t] and 7 are constants and 5* is the process for the unaffected stock price level. 
We can identify the following impact terms: 

• T] x(t) is the instantaneous impact of trading dx(t) = x(t)dt shares in [t,t + dt) 
and only affects the [t, t + dt) order. 

• 7(x(t) — Xq) represents the permanent impact price that has been accumulating 
over [0, t\ by all transactions up to t. 

2.1 ABM and GBM dynamics for the unaffected stock S 

The papers P and [2] assume that the unaffected (forward) stock price process S 
follows an ABM: 

dSt = aSQdWt, So (ABM) 

where W^ is a standard Brownian motion. 
Gatheral and Schied [3| assume instead that 

dSt = aStdWt, So (GBM). 

The absolute volatility in the ABM model is supposed to be proportional to Sq 
for dimensionality issues and for comparison with GBM. However the key fact is 
that the absolute volatility in ABM is constant and equal to a = aSo whereas in 
GBM the absolute volatility aSt is level proportional. 

One can easily compute the probability P that the asset is negative. One obtaines 

P(5i < 0) = $(-(aVt)-i) (ABM), F{St < 0) = (GBM) 

Clearly, ABM has the problem, especially for long times, that the unaffected 
stock probability of going negative is large. In the limit for large t such probability 
becomes 1/2. Luckily, trade execution is usually considered on short time scales, 
so that this problem is not as bad. However, GBM is a safe improvement since it 
immediately excludes negative stock values completely. 

2.2 DD dynamics for the unaffected stock S 

Consider now the DD model. We formulate such model for the unaffected stock 
price S as follows: given a GBM Yt with volatility parameter a, 

dYt = aYtdWt, Yo = So-K 
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we define S as 

St = K + Yt or equivalently dSt = a{St - K)dWt, So (DD) 

In other terms, a displaced diffusion model for S means that S* is a shifted GBM. 
The shift is a constant K and the GBM is Yf. To keep the asset positive, we clearly 
would need to assume K > 0. 

If i^ < 0, the probability of a negative stock S can be computed as 

ns. < 0) = * (zMl^^Z|)±^f!Z!) (dd w.th k < o) 

One can clearly see that the probability of having a negative stock in the long 
run, when the shift is negative, tends to 1, so that for negative shifts the asymptotic 
probability of negative states is worse than for (ABM). 

However, in general the DD model has the advantage of mimicking the key 
features in the GBM and ABM in a single model. Indeed, consider the absolute 
volatility in the DD model: 

a{St-K) = aSt{l-K/St) 

For large stock values S, the volatility looks like the level-proportional absolute 
volatility in a GBM model. Indeed, for large S the ratio K/S becomes small and 
one has 

aStil - K/St) ^ aSt for St » K. 

Instead, for small values of S, the absolute volatility looks more similar to the 
constant absolute volatility of the ABM, in that 

a{St -K)^ <y{-K) for St ^ 0. 

This is usually summarized by practioners by saying that "DD behaves like ABM 
for small stock values and as GBM for large ones" . 

2.3 The trade execution problem 

We consider a trade execution strategy, described by the above absolutely continuous 
stochastic process 

1 1-)- x{t) 
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representing the amount of shares that are left to be sold at time T. We assume 
that the process is adapted. Our aim is completing the sale of X shares by time T, 
so that our boundary conditions are 

a;(0) = X, x{T) = 0. 

The trade execution problem consits in finding the x as above that minimizes an 
objective function based on costs and risk terms. We look at those terms now. 

2.4 Cost Function 

The instantaneous cost of the strategy is given by the cost of buying dx{t) = x(t)dt 
shares at time [t,t + dt) at the impacted stock price St, spending Stx(t)dt. 

We can then write the costs arising from the strategy, following Gatheral and 
Schied [3] as 

C{x):= I Stx{t)dt= I [St + r]x{t)+^{x{t)-xo)]x{t)dt 
Jo Jo 

T rT 

2j^ , ,.v2 



= -XSo- x{t)dSt + v {x{t)ydt + -fXy2 
Jo Jo 

where we have used an integration by parts. This calculation does not involve 
the specific dynamics of S yet and is general. 

2.5 Risk Function: Value at Risk 

A common measure of risk is Variance. However, Gatheral and Schied [3] argue that 
for GBM a different interesting choice is a Value at Risk (VaR) risk based function. 
VaR is basically a given quantile of the difference between the value of the position 
at the time of calculation and the future value of the position at the risk horizon. It 
is a quantile of the future loss associated with the position. Let //„ ^ /^ be the Value 
at Risk measure computed at time t, for the position, for a given confidence level a 
(for example a = 0.95 or a = 0.99) and oven a time horizon h. We may omit some 
arguments when clear from the context. This means that by definition, given the 
market filtration 5" and given the probability measure P, we have 

F{St - St+h < i^a,t,h\3't} = a- 
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If at time t we have x{t) shares with price St, the time t VaR measure for a risk 
horizon h for that x{t) position under the DD dynamics at confidence level a would 
be 

ut[x{t){St - St+h)] = x{t)ut[{St - St+h)] = x{t)ut[{Yt + K- Yt+n - K)] 

= x{t)vt[{Yt - Yt+h)] = x{t)vt[Yt{l - exp(-a2/i/2 + a{Wt+h - Wt)))] 

= x{t)YtUt[l - exp{-a'^h/2 + aVhe)] = x{t)Ytqa[l - exp{-a'^h/2 + aVhe)] 

where e is a standard normal random variable, where we have used the homogeneity 
property of VaR, and where in general qa{X) is the a quantile of the distribution of 
the random variable X. Define the constant 

Aa := qa[^ — exp(— (j^/i/2 + aVhe)] 

to obtain that 

iyt[xit){St - St+h)] = Mt)Yt = Xx{t)iSt - K). 

This is the VaR measure for the instantaneous position at time t. If we average VaR 
over the life of the strategy we obtain the risk criterion 

i?VaR„(^) ._l x{t){St - K)dt. 
Jo 

This will be one of our measures of risk, to be included in the criterion to be 
minimized for optimality of the execution strategy. 

2.6 Risk Function: Expected Shortfall 

Expected shortfall iJia,t,h at time t for the horizon t + /i is defined as the expected 
value of the loss of the position beyond the VaR level, after the given time horizon, 
and for a given confidence level (the one used for VaR, say a). Hence it is defined 
as 

l^a,t,h '■= EjS'i — St+h\St — St+h > ^a,t,h, '3^t\- 

Expected Shortfall is usually preferred to VaR on the basis of the fact that it 
inspects the whole tail beyond the chosen quantile. Moreover, expected shortfall is 
a coherent risk measure, namely it is sub-additive, while Value at Risk is not. 
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By proceeding analogously with the VaR case, if at time t we have x{t) shares 
with price St, the time t expected shortfall measure for a risk horizon h for that x{t) 
position under the DD dynamics at confidence level a would be 

fit[x(t){St - St+h)] = x{t)Yt ma[l - exp{-cr^h/2 + aVhe)] 

where e is a standard normal random variable, where we have used the homogeneity 
property of expected shortfall, and where in general ma{X) is the expected shortfall 
for a loss given by the random variable X beyond the a quantile of the distribution 
of X. Through straightforward calculations one obtains 



1 



TTlr 



1 — a 



$ 



ln(l - K) + (y^h/2 



a 



Vh 



$ 



ln(l - K) - (r^h/2 



a 



v^ 



=:A. 



and hence 



flt[xit){St - St+h)] = x{t)Yt \a = \aXit){St - K). 



This is the expected shortfall measure for the instantaneous position at time t. 
If we average this measure over the life of the strategy we obtain the risk criterion 

R^^"{x) := A / x{t){St - K)dt. 
Jo 

This will be one of our measures of risk, to be included in the criterion to be 
minimized for optimality of the execution strategy. 



2.7 Cost plus Risk 

We are now ready to define our criterion to be minimized. We put together costs 
and risks in a single criterion 



E [C{x) + R{x)] = E 



C{x) + L\ I x{t){St-K)dt 





In this criterion: 



L is a cost/risk leverage parameter that measures risk aversion in executing the 
order. The larger L, the more important the risk component of the criterion 
compared to the cost component; 
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A is either A or A depending on whether we are using the VaR or expected 
shortfall risk function, leading to 



E [C{x) + R{x)] = -XSo+-fXy2+riE 



/ xitfdt + L- I xitMSt - K)dt 
Jo V Jo 



= -XSo + 7^72 + V^ 

where we have set A = X/rj. 
The problem is now finding 



T pT 

2, 



xitydt + LX / x{t)Ytdt 
Jo 



X* = arginf^.E 



T rT 

2, 



xitydt + LX I x{t)Ytdt 
Jo 



(2.1) 



3 Solution of the trade execution problem for a 
displaced diffusion 



The above problem ( |2.1 ) has been solved by Gatheral and Schied |3]. If we pretend 
for a moment that Y is our true unaffected underlying stock, the criterion we have is 
the same as the criterion for a GBM Y, and this has been solved in [3j. Indeed, The- 
orem 1 in Gatheral and Schied [3] implies that the unique optimal trade execution 



strategy for the problem (2.1) above is 

. T -t ' 






T 



X-LX- I Y,du 







with the minimized value being given by 

E 



T rT 



x*{tydt + LX / x*{t)Ytdt 
uo Jo 



X2 XLTXYo {LXf 



v2 ( g.^r _ 1 _ ^2y _ 



4t^2 



a^T- 



T ' 2 8^6 

which, replacing Y hj S — K, leads to the following 
Theorem 3.1. (^Optimal execution strategy for a displaced diffusion). The 



unique optimal trade execution strategy attaining the infimum in (2.1) is 

. T-t 



X, 



T 



X - XL J I {S^ - K)du 
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Furthermore, the value of the minimization problem in (2.1) is given by 



E 



f x*{tfdt + \L I x*{t){St-K)dt 
Jo Jo 



X-' ^ LXTXiS,-K) (^^^^ ^^,j, / ,,, _ J _ ^,^ a'T^ 



T 2 8a^ ' ' \ 2 

4 DD under VaR or ES escapes Schied's invari- 
ance result 

A quite interesting invariance result due to Alexander Schied is proven in (H]. The 
result is rather general, and here we are concerned with a special case. The paper 
[6] considers the risk criterion 

R^^{x) := A / x{t)Stdt 
Jo 

which would be a VaR or ES type risk function when the unaffected price S follows 
a geometric brownian motion (as proven in [3]). 

Then the paper [6] shows that as long as S follows a martingale, the optimal 
trade execution strategy x* does not depend on the specific martingale dynamics 
one chooses. Since the displaced diffusion model we consider in this paper is also a 
martingale, one does then deduce that under the R^^ risk criterion and DD dynamics 
we obtain the same solution we would obtain under a GBM dynamics. 

This seems to contradict our result in this paper, given that our solution is 
different from the GBM one (as we have K in our solution). In fact, however, there 
is no contradiction. 

Indeed, we use a different risk criterion, since we use 

A / x(t){St — K)dt instead of A / x(t)Stdt. 
Jo Jo 

In the specific case of DD, the first criterion corresponds to the VaR or Expected 
Shortfall (depending on the choice of A) of the position, whereas R^^ corresponds 
neither to VaR nor to ES. Hence if we include K in our criterion we are using a risk 
criterion that is different from |6] and this is why we escape the invariance result 
and our optimal solution is different from the GBM one. 

We believe it is good to retain K in the criterion for the DD case because, in 
this way, we are really using VaR or expected shortfall as a risk criterion. 



D. Brigo, G. Di Graziano, Optimal Execution under Displaced Diffusion and Comparisons 11 

5 Optimal solution for displaced diffusion processes 
under an alternative square risk function 

In order to test the robustness of the optimal strategy, we introduce the alternative 
" squared-asset expectation" (SAE) risk criteria 

R^'^^ix) = X [ x\t)a^E[S^]dt. (5.1) 

Jo 

This goes a little back towards considering the variance as a measure of risk. 

Remark 5.1. (On the meaning of A). It is important to notice that while in the 
VaR and expected shortfall A had a precise value emerging endogenously from the 
chosen criterion, here we have just introduced A exogenously, as a parameter we are 
free to set in our criterion. 

The criterion above is similar to the one used by Almgren and Chriss in [2] but 
with the displaced diffusion replacing the arithmetic Brownian motion (ABM) as 
far as the dynamics of St are concerned and the expected value of the square of the 
price process replacing the square of the price process itself. Note that the term 

9{t) ^ E[S^] 

is increasing in t, whereas the corresponding term in |2] is constant and equal to 
Sq. Moreover in the Almgren and Chriss framework the volatility of the asset a is 
expressed in units of the reference currency, whereas in this paper a is a percentage 
number. Using the notation of this paper, the risk criteria employed in p] can be 
re-written as 

R^^ = X x\t)a^S^dt. (5.2) 



Everything else being equal, the magnitude of the risk component (5.1 ) is higher as 



time goes by compared to (5.2). Finally note that the expected cost is the same 
whether St is modelled by a DD or an ABM, so the optimal solution depends only 
of the choice of risk function. 

The optimal execution strategy can be derived by solving the optimisation prob- 
lem 



X* = arginf^jE 



T fT 



x{tydt + LX I a^x\t)g{t)dt 





(5.3) 



D. Brigo, G. Di Graziano, Optimal Execution under Displaced Diffusion and Comparisons 12 

Using calculus of variation, we show in appendix |X] that the optimal solution needs 
to satisfy the ODE 

x(t) = egit)xit), (5.4) 

where 

k = avLX, 

and the initial and terminal conditions are given by x(0) = X and x{T) = re- 
spectively. A solution to the boundary value problem above could be found using 
standard numerical routines. However, in appendix |B] we provide a simple and 
efficient alternative algorithm based on series expansions. 

6 Comparing Risk Criteria 

If we are to compare optimal strategies under the SAE risk criterion with the VaR 
and expected shortfall criteria strategies, we need to set A in the SAE criterion in 
a way that makes the comparison sensible. We call this procedure "equalizing the 
A's". 

A first possible way to equalize A's is to match the initial and final values of the 
instantaneous criteria before integration. We set 

A xmSo-K) = X'AiM,DD x\0) a'Sl A x{T) (St-K) = X'//m,dd x'{T) a^Eo[5^] 

for the cases ABM and DD respectively (as per A's index). 

Recalling that a;(0) = X and x{T) = 0, we have that the second equations (for 
T) are just identities, whereas the first ones lead to 

\SAE _ \ Sq — K 

so that if we aim at comparing the SAE and the VaR or ES criteria outputs we may 
use the above A in the risk criterion for SAE. 

A second possible way to derive A in SAE is to check what happens for the 
VWAP solution of the simpler case 

T — t 
Xo{t) = X^^ 

and match the risk functions corresponding to this solution. 
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Straightforward calculations based on solving 



Eo 



T 



A xo{t) {St - K)dt 



Eo 






in the DD case, while using the expansion e'^ ~ 1 + at, leads to 

\SAE _ _J^ I /^ . N 

{So - K) (^^ + ^^So-K)^) 

This is the formula we will use, and we will set also 

\SAE ._ ^SAE ((- 9N 

6.1 Comparative Results 

In this section we present a series of numerical results which compare the different 
models and risk criteria analysed in the paper. In order to stress the main features 
of each model we chose, in some instance, parameters which may be regarded as 
extreme or unrealistic. We shall look at optimal solution under a more realistic set 
of parameters, and in particular the execution horizon, later in the section. 

In the graphs below, we compare the optimal solution of the following dynam- 
ics/risk criterion combinations: ABM+SAE, DD+SAE, GBM+VaR and DD+VaR. 



The A's will be equalized through Equations (6.1) and (6.2). 

In figure [ij the time horizon for the execution is set to one year, which is ad- 
mittedly very long. The volatility of the asset is 20% and the DD shift K is 20% of 
the initial asset price. Most importantly the cost/risk parameter L is relatively low 
at 100 and the temporary market impact parameter t] has been set to 1%. For the 
choice of parameters above, the optimal execution strategy is very similar in all the 
models and almost linear in t. 

Figure [2] is based on the same parameters of figure [l] with the exception K which 
is equal to 80% of the initial asset price. Note that the instantaneous volatility of 



the DD is rescaled using formula (6.3) to ensure that the integrated volatility of the 



DD and GBM models are roughly of the same order of magnitude. 



no _ nBM ^0 



^UU^^GBM_^_ (6.3) 

Since the paths generated under the GBM and DD are different, the two optimal 
solutions differ, as we show with a GBM path which having a stronger performance 
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than DD. The optimal execution pohcy under ABM and DD using the SAE risk 
function is not path dependent and it is virtually unchanged under the new param- 
eter choice. 

Figure |3] shows that under a combination of high volatility and high L, the 
optimal strategy under the VaR risk criteria can change sign, e.g. a selling policy 
may turn into a buying policy before T. This is due to the integral of Su or Yu in the 
optimal solution. When this term is relatively big compared to X then the amount 
of asset to execute x{t) turns negative. Mathematically this is due to the lack of a 
smoothness constraint at for the process x(t). 

Figure |4] shows how an increase in the market parameter t] slows down the 
execution policy. 

In figure |5| we reduced the time horizon to a more realistic 1 day. Even if we 
choose very high risk parameters, e.g. L = 5000, a = 99% and rj = 0.1, the optimal 
solution for all models is very close to linear. 

Finally, note that the solution under the VaR risk function depends on the rel- 
ative magnitude of X and the price of the asset. In order to illustrate this feature, 
in figure |6] we use the parameters of figure [3] but we increase the units to execute 
to 10000 (i.e. 100 times the initial asset price). Notice how the optimal solution 
becomes approximately linear regardless of the path. The solutions obtained using 
the SAE risk function are also linear because in the rescaling formula for the rele- 
vant A parameter, the X appears in the denominator, making the risk contribution 
negligible. 

Since in practical applications the execution horizon is rarely longer than a few 
days (more commonly it is measured in hours), the VWAP execution strategy will be 
a very close approximation of the optimal strategy obtained using more sophisticated 
models. 

Remark 6.1. (Actual implementation and estimation of the asset dynam- 
ics). In real world applications we observe only one realised path of the "true" asset 
dynamics and the solution of the two GBM/DD models will depend on the "cali- 
brated" parameters a^^'^^ / {a^^,K). Since the path on which both models are to 
be estimated will be the same, the optimal solution of the two models will be much 
closer than in the example above where the paths differ (note however that the same 
brownian shock process W realization has been used to sample both DD and GBM). 
If we had to be really precise, we would have to simulate a path using a third, more 
realistic model, and then calibrate the GBM and DD parameters to such path. At 
that point we would apply the optimal solution formula for the 2 models and com- 
pare the two solutions. Although the calibration of the GBM to historical data is 
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quite simple, once the drift is assumed to be zero, calibrating the DD model (and 
K in particular) is slightly more complicated, and may account for higher horder 
moments. 

This method would be more ideal, but we may anticipate that results would be 
even closer than in our examples above, showing probably even smaller departures 
from the trivial VWAP solution. 

An alternative calibration to fit the parameters could be an option pricing cal- 
ibration, obtaining the dynamics statistics under the pricing measure, rather than 
the historical one, but for the optimal trade execution problem we need to be careful 
about mixing historical measure and risk neutral measure statistics. 



7 Conclusions 

In this paper we solved the optimal trade execution problem under the Value-at-risk 
or Expected shortfall based criteria of Gatheral and Schied when the underlying 
unaffected stock price follows a displaced diffusion model, encompassing arithmetic 
and geometric Brownian motion type features at the same time. We further intro- 
duced an alternative risk criterion called Squared-Asset Expectation, and we noticed 
that the optimal trade execution strategy does not depend strongly on the specific 
risk criterion one adopts. We find that the optimal solution is close to the basic 
linear Volume Weighted Average Price (VWAP) solution regardless of the specific 
diffusion dynamics or risk criterion that is chosen, especially on realistic trading 
horizons of a few days or hours, and unless unrealistically large values of the risk 
or dynamics parameters are taken. Our preliminary findings suggest the conjecture 
that more general dynamics, beyond simple diffusions, need to be considered to see 
a relevant impact on the optimal strategy. A similar conclusion holds also for the 
risk criterion. 
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Appendices 

A Derivation of the ODE for the Optimal Strat- 
egy 



The optimisation problem (5.3) requires that we find the path x{t) which minimises 
the loss function 



T pT 

2j. , r \ / _2„2/ 



(A.l) 



iJo Jo 

In order to do so, consider the following perturbation of the processes x{t) and x{t) 



x'{t) = x{t) + eh{t) 
x\t) = x{t) + eh{t) 

where h{t) is an arbitrary function satisfying h{0) = h(T) = and e is a real 



constant. Substituting the perturbed path into (A.l) we obtain 



Hie) 



E 



T 



[x\t)Ydt + L\ a\x'{t)Yg{t)dt . 
IJo Jo 

Taking the first derivative with respect to e and evaluating the resulting expression 

at e = 0, it follows that 



if (0) = E 



2 / x{t)h{t)dt + 2L\ j a^x{t)h{t)g{t)dt 
Jo 



(A.2) 



Using the integration by parts formula and the constraint h{0) = h{T) = 0, we 
obtain the equality 



x{t)h{t)dt 







x{t)h{t)dt, 



which allows us to simplify (A.2) as follows 



H{0) = E 



{x{t) - LXa^x{t)g{t)) h(t)dt 



The optimal path x{t) is obtained by setting H{Q) = 0. Since the function h{t) 
was chosen arbitrarily, the following must hold for all < t < T 

x{t) - L\a'^g{t)x{t) = 0. 
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B Numerical algorithm for optimal solutions 

Consider the ODE 

x{t) = k^g{t)x{t). 

Simple manipulations show that g{t) can be written as follows 



where 
and 



git) 



y2 
^0 


( — ] + 2 — + expiah) 


(B.l) 


PK + e'^'*) 


(B.2) 



P = {Yokf 



ao 



K 



K 

Assume that the solution of the ODE can be expanded in a Taylor series around 
t = 



X 



j=0 



(B.3) 



The coefficient of the series can be found by substituting (B.3) into the ODE, 



Y,U + 2)(j + l)aj+2t' = P (tto + exp{vt)) ^ ajt^ . 

j=0 \j=0 J 

Substituting the exponential term with its Taylor expansion in t and exchanging 
the order of summation we obtain 



vjr=0 / i=l j=0 ^' 



oo oo 



0=0 / j=l j=i 
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Truncating both sums at m and equating the coefficient of the same order in t 
on both sides of the equation we obtain a recursive expression for the coefficients aj, 

"j+2 = (_^- + 2)(j + l) ( ^^ "^ "°''"^' ^ ^ ^^^''kl 

for j = to m — 2. 

The expression above can be written in a more compact form using vector nota- 
tion 

where we have defined 

Aj = (ao,ai,. . . ,aj) 

T/ - f^ ^^ ^^ 

Using the initial condition a;(0) = X it follows that oq = ^- On the other hand, 
Ci can be derived numerically via a simple root search using the terminal condition 

x{T) = 0. 
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Figure 1: Model Comparison: So = 100, X = 100, T = 1, K = 20, a = 0.2, r] = 1, 
L = 100 
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Figure 2: Model Comparison: So = 100, X = 100, T = 1, K = 80, a = 0.2, r] = 1, 
L = 100 
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Figure 3: Model Comparison: So = 100, X = 100, T = 1, K = 20, a = 0.75, r] = 1, 
L = 500 
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Figure 4: Model Comparison: So = 100, X = 100, T = 1, K = 20, a = 0.75, r? = 5, 
L = 500 
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Figure 5: Model Comparison: 5*0 = 100, X = 100, T = 1/250, K = 90, a = 0.99, 
rj = 0.1, L = 5000 
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Figure 6: Model Comparison: So = 100, X = 10000, T = 1, K = 20, a = 0.75, 
r^ = 1, L = 500 
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